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Abstract. We show that there exist solutions of drift-diffusion equations in two dimensions with 
divergence-free super-critical drifts, that become discontinuous in finite time. We consider classical 
as well as fractional diffusion. However, in the case of classical diffusion and time-independent 
drifts we prove that solutions satisfy a modulus of continuity depending only on the local L 1 norm 
of the drift, which is a super-critical quantity. 



1. Introduction 

We study the continuity of solutions to equations with divergence-free drift and fractional or 
classical diffusion. We prove that in supercritical regimes, there are solutions which become dis- 
continuous in finite time. However, we also prove that in two dimensions, a solution to a drift 
diffusion equation with classical diffusion stays continuous if the drift is only locally bounded in 
L 1 (which has supercritical scaling), provided the drift is time-independent. 

Equations with drift and diffusion appear in numerous places in mathematical physics. In many 
cases, the drift depends on the solution and the equation is nonlinear. A successful understanding 
of well-posedness of the problem in each case depends on the a priori estimates that can be estab- 
lished. In most cases, these are based on the linearized drift-diffusion equation, which provides the 
motivation for this work. 

We consider the problem of continuity of solutions to the Cauchy problem 

d t 6 + u- W+(-A) s = O (1.1) 
0(0,0 = 00 (1-2) 

where s G (0,1], and u is a given divergence-free vector field. 

For each s > 0, the equation has a natural scaling: if 8(t,x) is a solution of Ql.l| ) with drift 
u(t,x), then 9\(t,x) = \ 2s ~ 1 9(\ 2s t, Ax) is a also a solution of ( |1.1[ ), but with drift given by 
u\(t, x) = \ 2s ~ 1 u(\ 2s t, Ax). A Banach space X, with norm || • is called critical with respect to 
the natural scaling, if ||w\||x = \\ u \\x for all A > 0. If on the other hand, ||«a|U — ^ oo as A — )• 0, 
the space is called supercritical. In the supercritical cases, when one zooms in at a point (i.e., 
sends A — >■ 0), the bound on the drift becomes worse, so that regularity of the solutions cannot be 
inferred from linear perturbation theory. In view of the scaling described above, for s E (0, 1/2) a 
critical space for ( |1.1[ ) is the Holder space C l ~ 2s , for s = 1/2 it is the Lebesgue space L°°, while 
fors G (1/2,1] it is L d ^ 2s ~ v > . 

In the context of fluid mechanics, the case of divergence-free drifts is of special importance due 
to incompressibility, while (fractional) diffusion is a regularizing term, for instance, in the well 
known surface quasi-geostrophic (SQG) model lfT3l [T51 . The possibility of finite time blowup 
for the SQG equation with supercritical fractional diffusion is an outstanding open problem. One 
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could speculate that the divergence-free character of the drift plays an important role in the well- 
posedness of the problem. In fact, blow up in finite time does not seem to be known to hold not just 
for SQG, but for most (if not all) of the supercritical active scalar equations with divergence-free 
drift currently in the literature (see, for example [|71[8l[T0l|22]l). In contrast, for some drift-diffusion 
equations with non-divergence-free drifts, a blow up scenario is possible, and well understood. 
For example, this is the case in Burgers equation with fractional diffusion 03l26ll, the Keller-Segel 
model [EH, and many other equations |fl2llT7ll29ll . Even self-similar blow up may be sometimes 
obtained, since there is no mechanism which prevents mass-concentration. 

Indeed, the divergence-free condition on the drift is known to imply some qualitative properties 
of the solution. For pure transport equations without diffusion, the flow is well defined almost ev- 
erywhere just assuming that the drift is in the Sobolev space W 1 ' 1 instead of the classical Lipschitz 
assumption of Picard's theorem (see [|2l and Il20]| ). Also, certain type of singularities are ruled out 
for divergence-free drifts, as in [fT8ll . We give another example of this phenomenon in Section |5} 
where we prove that, in two dimensions, non-vanishing continuous divergence-free vector fields 
(not necessarily Lipschitz or even Holder) have unique trajectories (see also [4]). For equations 
with drifts and diffusion, the divergence-free condition has been used to obtain some estimates 
which are independent of the size of the drift, e.g. first eigenvalue estimates (3), mixing rates ifTTll . 
and expected exit times [|24l . 

For scalar equations with drift and classical diffusion (s = 1 in ( |1.1[ )), if the drift u is assumed to 
be divergence-free and in the critical space L°°(BMO~ l ), then one can obtain a Holder estimate 
on the solution 9 by extending the methods of De Giorgi, Nash, and Moser (see ||23l or ||3TTl ). If 
the drift were not assumed to be divergence-free, then one would obtain Holder continuity of the 
solution under the stronger assumption u E L°°(L n ). Note that the space L^^BMO -1 ) is larger 
than L°° (L n ), but it has the same scaling. Therefore, the divergence-free assumption only provides 
a borderline improvement in this result. 

In their celebrated paper [6], Caffarelli and Vasseur were able to prove well-posedness of the 
critical SQG equation based on an a priori estimate in Holder spaces for ( |1.1[ ) when s = 1/2 and u 
is divergence-free and in L°°(BMO) (well-posedness of SQG was also proved independently by 
Kiselev, Nazarov, and Volberg Il27l0 . Other proofs of this result were given in [|25l and lfT4ll . For 
non-divergence-free drifts, the same type of Holder estimate was obtained in ||33ll by a different 
method assuming u E L 00 (L 00 ). Again, the divergence-free assumption only provides a borderline 
improvement in the regularity result since BMO and L°° have the same scaling properties. 

In lfT6ll . Constantin and Wu investigated lower powers of the Laplacian in the diffusion using 
techniques from [6J and [5]. They obtained a priori estimates in Holder spaces for the equation 
< |1.1| ) where s E (0, 1/2), u is divergence-free and in L°°(C 1_2s ). Using the ideas from 11331 . 
the result was generalized to non divergence-free drifts in [|32ll . These estimates do not suffice to 
show well-posedness of the surface quasi-geostrophic equation in the supercritical regime. One 
might wonder whether the result in |[T6l could be improved using the divergence-free condition in 
a stronger way. In fact, in [J9]|, it was suggested that the solution of for any s E (0, 1/2) would 
be Holder continuous just assuming that u E L°°(BMO) and is divergence-free. We disprove 
this last statement here. We show that the result in [16] is in fact sharp by proving that for any 
a < 1 — 2s, there is a divergence-free drift u E L°°(C a ) for which the solution of ( |1.1[ ) develops 
a discontinuity starting from smooth initial data. 

We now state our main results. In this paper a modulus of continuity will be any continuous 
non-decreasing p : [0, oo) — > [0, <x>) such that p(0) = 0, and we say that 9 : (t ,tx) x W l — > R 
breaks (satisfies) p at time t if there are (no) x,y E W 1 such that \9(t,x) — 9(t,y)\ > p(\x — y\). 
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Although we will restrict our considerations to the case of two dimensions d = 2 here, several of 
our results extend to more dimensions. 

Our first result shows that the result of Constantin and Wu IfToll is sharp, even for time-independent 
drifts. 

Theorem 1.1 (Case s < 1/2). Let s G (0, 1/2) and a 6 (0,1- 2s). There exist a positive time 
T and a smooth function 6q with ||#o||c 2 (r 2 ) — 1 such that for any modulus of continuity p, there 
exists a smooth divergence -free time -independent vector field u with ||tt||o(K2) < 1 such that the 
smooth solution of ( ] 1 . 1 1 ) — ( |TT2| ) breaks the modulus p before time T. 



For the case s > 1/2, a critical assumption on the drift can be given in terms of Lebesgue spaces. 
It is conceivable that the method of Caffarelli and Vasseur [6| can be extended to s > 1/2 in 
dimension d assuming that u G L 00 (L d ^ 2s ~ 1 )), although this has not been written down anywhere, 
to the best of our knowledge. We prove that no weaker assumption on u could assure continuity of 
the solution. 

Theorem 1.2 (Case 1/2 < s < 1). Let s G [1/2, 1) and p G [1, 2/(2s — 1)). There exist a positive 
time T and a smooth function 6 with ||#o||c 2 (r 2 ) — 1 such that for any modulus of continuity p, 
there exists a smooth divergence -free time -independent vector field u with \\u\\lp(^?^ < 1 such that 
the smooth solution of ( ] 1 . 1 1 ) — ( |TT2] ) breaks the modulus p before time T. 



The results obtained in Theorems [TT| and 1 .2 for smooth drifts can be used to prove that there are 



divergence-free time-independent u G C a (with a < 1 — 2s, if s G (0, 1/2)), respectively u G L p 
(with p < 2/(2s — 1), if s G [1/2, 1)), with distributional solutions of the initial value problem 
( |1 . 1| ) — ( |TT2| ) evolving from smooth initial data, which fail to be continuous at any t > 0. Indeed, the 
drifts for different ps are constructed from a p-independent non-smooth drift using p-dependent 
cutoffs near the origin (where p will be broken) to ensure smoothness. Removing this cutoff at the 
origin will result in a (limiting as cutoff area shrinks to the origin) distributional solution which 
breaks any modulus p in finite time. Moreover, this is true for any time t > thanks to infinite 
speed of propagation of diffusion. We believe that by following the ideas in [El |20j , one can show 
that these distributional solutions are unique whenever the divergence free drift lies in L 1 (BV) 
(and all the drifts considered in our Theorems l.lf 1.3 have this regularity). 



One might ask what happens in the endpoint case of classical (and local) diffusion s = 1, and 
the answer is quite intriguing. First, we show that if we allow the drift to be time-dependent, then 



the above results continue to hold. (We note that the remark after Theorem [L2] also remains valid 
in this case, this time after the removal of the temporal cutoff near the "blow-up" time t q > from 
the proof, albeit with breaking of all moduli guaranteed only by time t q .) 

Theorem 1.3 (Case s — 1: time-dependent drifts). Let s = 1 and p G [1,2). There exist 
a positive time T and a smooth function 9q with ||^o||c ,2 (r 2 ) — 1 such that for any modulus of 
continuity p, there exists a smooth divergence-free vector field u with \\ sup t |u|||lp(r 2 ) < 1 such 
that the smooth solution of ( |1.1[ )-( [L~2"| ) breaks the modulus p before time T. 



In the case of time-independent drifts, however, Theorem 1.3 is surprisingly false! We prove 
that the solution 9 has a logarithmic modulus of continuity which depends on u via its local (super- 
critical) L 1 norm only and, in fact, continuous distributional solutions exist for non-smooth locally 
L 1 drifts. This is a remarkable property which holds in two space dimensions only. 

Theorem 1.4 (Case s — 1: time-independent drifts). Let s = 1 and assume that u G L 1 1 oc (M 2 ) 
with \\u\\ L i ( K 2) = sup^^ 1 1 « 1 1 l 1 ^^)) < oo is a divergence-free time-independent vector field. If 
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#o E C 2 (IR 2 ) fl W 4,1 (M. 2 ), then there is a distributional solution of < \ 1 . 1 | )-( [L2] ) which is continuous 
and at any time t > satisfies a modulus of continuity given by 

, , ^(i + IHIzj )||^ollcw,i(i + r 1 ) 

Pt(r) = loc r=, (1.3) 

for r E (0, 1/2), with some universal C > 0. 

We remark that if instead of u E L\ oc we assume u E L 1 , then one may lower the regularity 
assumption on the initial data to 9 E C 2 fl W 2 ' 1 (see the proof). 



We note that the last claim in Lemma [6T2] shows that for each t > 0, this solution satisfies the 
(spatio-temporal) modulus p to on (t , oo) x IR 2 as well. Moreover, the result in fact holds for any 
distributional solution which is a locally uniform limit of smooth solutions with drifts converging 
to u in L\ oc . In particular, if u is smooth, the (unique) solution of the Cauchy problem satisfies the 
modulus of continuity ( |1.3| ), which depends on the drift via the super-critical norm \\u\\ L \ only. 



An analogous result in the elliptic case was proved in [31 ]: the elliptic maximum principle plus 
an a priori estimate in H 1 , which hold for solutions of the PDE, suffice to show that a function has 
a logarithmic modulus of continuity. This idea can, in fact, be traced back to Lebesgue [|28l ! 

In the parabolic setting the situation is somewhat different. The parabolic maximum principle 
plus the energy estimates do not suffice to show the continuity of the solution. The following ex- 
ample illustrates the difficulty: if 6(t,x) = <p(x) forsomey? E H^R 2 ) \C(R 2 ), then 9 E C°°(H l ) 
and it satisfies the parabolic maximum principle, without being continuous. In order to overcome 
this difficulty, we need to use the equation to prove that for each time t, the elliptic maximum 



principle holds modulo an error that we can control (see Lemma 6.8). A crucial ingredient will 
also be that d t 9 solves the same equation as 9, which allows for some important bounds on this 
quantity (see Lemmas 6.2 and 6.3[ ). Of course, this only holds when u is time-independent. 



Inspired by |[T9l , we also explore a slightly supercritical equation. This is the case in which the 
fractional Laplacian is replaced by an integral kernel which is logarithmically supercritical. 

Theorem 1.5 (A slightly supercritical case). Let m : R + — > R + be a smooth non-increasing 
function such that 

r^dr<oo (1.4) 
Jo 1 + ^ 

and rm(r) is non-decreasing on (0, 1). There exist a positive time T and a smooth function 9 with 
||0o||c 2 (]R 2 ) < 1 suc h that for any modulus of continuity p, there exists a smooth divergence-free 
time-independent vector field u with ||w||z,o°(]R2) < 1 such that the smooth solution of 



d t 9 + u-V9 + P.V. [ ( 9(x) - 6(x + yj) dy = (1.5) 

v ' \y\ 

with initial condition 9q breaks the modulus p before time T. 



This result suggests that in order to hope for continuity of solutions to (1.5), one should not 



depart from the critical case m(r) = 1/r by "more than a logarithm". It would be interesting to 



show that for generic divergence-free L°° drifts, solutions to ( 1.5 ) are continuous if the integral in 



( |1.4[ ) diverges. In fact, for the dissipative Burgers equations it was shown in [19| that when ( |1.4| ) 
holds, shocks develop in finite time, while if the integral in ( |1.4[ ) diverges, global regularity holds. 

The article is organized as follows. In Section|2]we present informally the main idea behind the 
loss continuity in finite time of solutions to < | 1 . 1 [ ) (for simplicity we take u E L°° and s < 1/2). 
Section [3] contains the proofs of Theorems 1 1 . 1 1 and 1 1 .2[ while Theorem [Tj] is proven in Section|4] 



ON THE LOSS OF CONTINUITY FOR SUPER-CRITICAL DRIFT-DIFFUSION EQUATIONS 



5 



Sections [5] and [6] contain the proofs of the positive results in this paper, both for divergence-free 
time-independent drifts: uniqueness of particle trajectories for the transport equation with a con- 
tinuous drift (Theorem 5.1 ) and Theorem 1.4 We conclude by proving Theorem 1.5 in Section [7J 
LS was supported in part by NSF grants DMS-1001629 and DMS- 1065979, and an Alfred 
P. Sloan Research Fellowship. AZ was supported in part by NSF grants DMS-1056327, DMS- 
1113017, DMS-1 147523, and DMS-1 159133, and an Alfred P. Sloan Research Fellowship. VV 
was supported in part by NSF grant DMS-121 1828. 



2. LOSS OF REGULARITY FOR BOUNDED DRIFT AND SUPERCRITICAL DISSIPATION 

In this section we give a brief outline of our method of proving loss of continuity of solutions to 
drift-diffusion equations. In order to emphasize the main ideas, we present here the simplest case: 
the time-independent divergence-free drift is bounded and dissipation is super-critical with respect 
to the natural scaling of the equations (i.e., a = and s G (0, 1/2)). For the sake of simplicity 
of exposition, we will not require here the drift to be smooth and simply assume that we have a 
unique solution to ( | 1 . 1 1 ) — < p~72] ) . In the rigorous treatment in Section |3j we shall consider a smooth 
version of the drift defined below (and its generalization to other a) so that we need not worry 
about global existence and regularity of solutions. 

Define the stream function 



if>(xi,x 2 ) 



\x x - x 2 \ 



Fl + X 2\ 



and let u = V ip, where V = (— d X2 , d xi ). Then u e L°°(IR J ) is divergence-free (in the sense 
of distributions) and may be written explicitly as Let us consider ( |l.l| )-( fL"2| ) with this drift, s E 



U(Xi,X 2 ) 



' (0, —1) if x 2 > 

(0, 1) if x 2 < —\xi\, 

(1,0) if xi > \x 2 \, 

k (-l,0) if x 1 < -\x 2 \. 



\U II 
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i t t t\- 
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(0, 1/2), and some smooth initial datum 9 . 

Theorem 2.1 (Loss of continuity in the supercritical regime). Let s e (0, 1/2) and let p be any 

modulus of continuity. If 8 is non-negative in the upper half-plane, larger than 1 in Bi/ 2 (0, 1), 
and odd in x 2 , then the solution of < \ 1 . 1 [ )— ( fTT2| ) breaks the modulus p at some time £* < 1. 

Since p was arbitrary, it follows that solution loses continuity in finite time, which may be 
viewed as a blow-up in L°°((0, 1); C(R 2 )). 

Before proving Theorem 2.1 , let us observe two key properties of the system ( ] 1 . 1 h — ( p~T2] ) (again, 
these hold in the smooth case and we assume them to hold in the case at hand as well). 

Symmetry: If 9 is odd in x 2 , then so is the solution 0(t, •) for t > (because u{Rx) = Ru(x), 
where R(xi,x 2 ) = (x\, —x 2 )). 
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Comparison principle: If 9 is odd in x 2 and a super-solution of on the upper half-plane 
(and so a sub-solution on the lower half-plane) with sga(x 2 )9(0, x) > 0, then for t > we have 

sgn(x 2 )9(t, x) > 0. 

With these in mind, the intuition behind the proof of Theorem 27T] is as follows. Instead of 
( jl.l| )-( fL2"l ), consider a pure transport equation with drift u, no dissipation, and odd-in-x 2 initial 
condition satisfying 9 > Xb 1/2 (o.i) on the upper half -plane. Then the function 

9(t, •) = XB (x _ t)/a (0,l-t) - XB ( i_ t)/2 (Q,t-l) 

is a sub-solution, because the radius of the two discs forming its support decreases at the same 
rate as they approache the origin, so that the support stays in the set \x 2 \ > \xi\. So if 9 is the 
actual solution, then 9 — 9 is an odd-in-x 2 super-solution on the upper half-plane. The comparison 
principle (which also holds for the transport PDE) now shows that the oscillation of 9 on B X -t (0, 0), 
is at least 2, so that any modulus is broken before time t = 1. Adding now dissipation will decrease 
the supremum of 9(t, ■) on B^_ t y 2 (0, 1 — £). However, as the proof below shows, the latter will 
stay bounded away from zero on any finite time interval as long as the equation is supercritical. 

Proof of Theorem |2T7] Let 77 E C°°(R) be an even function, supported on [—1, 1], positive on 
(— 1, 1), non-increasing on IR + , and with 77(0) = 1. For some r E (0, l/\/2), we let 

4>{xx,x 2 ) = T)(\(xi,x 2 - 1) Ir" 1 ) - r)(\(xx,x 2 + 1) (r" 1 ), 

so that <p consists of a positive smooth bump centered at (0, 1) and with radius r, and a similar 
negative bump centered at (0, —1). Then for any s E (0, 1/2) there is c r>s > such that 

(-A) s (j)(x 1 ,X2) < c T , s (j)(xi,x 2 ) (2.1) 



holds in the upper half-plane {x 2 > 0} (the proof of this is given in Lemma 3.4 below). 
We now let 9 = (ft and 

9(t, x u x 2 ) = exp f-c r , s z(t)- 2s dr) fe, (2.2) 

for t E [0, 1), where 

z(t) = 1 - t. (2.3) 

Notice that z is the solution of the ODE 

z(t) =u 2 (0,z(t)) = -1, z{0) = l, 

which is the position of the original center (0, 1) of the positive bump, transported by the drift u. 
It is clear that 9 > for x 2 > and is odd in x 2 . The support of 9(t, ■) consists of two discs whose 
centers (0, ±z(t)) are transported towards the origin with the drift u, and have radii rz(t). That is, 
the support shrinks at the same rate it approaches the origin (and lies in the set \x 2 \ > \xi \ because 

r < 1/V2). 

We now claim that 9 is a sub-solution of ( |1.1[ ) in {x 2 > 0}. By scaling and ( |2.1[ ) we have 

(-A) s 9(t, ■) < c rjS z(t)~ 2s 9(t, ■) 
in {x 2 > 0}, so it suffices to prove 

d t 9 + u- V9 + c r , s z(t)- 2s 9 < (2.4) 
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in {x 2 > 0}. From the definition of 9, it is clear that we just need to verify 

d t <f) + + u ■ V0+ < 0, where <f> + (t, x ± , x 2 ) = r) ( K^flZ^M! 

\ rz(t) 

The point now is that <p + is transported by the vector field v(t, x) = —x/z(t), that is, it satisfies 

d t <p + + v ■ V0+ = 0, 

as one can see by a simple computation. We therefore only need to show (v — u) ■ V</> + > 
which, due to (fi + being supported in x 2 > \xi\, radially symmetric, and radially non-increasing 
with respect to center (0, z(t)), is equivalent to 

< ({x 1 ,x 2 )z(t)- 1 + (0,-1)) • ( Xl ,x 2 - z{t)) = z(t)- 1 \(x 1 ,x 2 - z{t))\ 2 (2.5) 

on the support of </>+. This clearly holds, so 9 is a sub-solution of ( |1.1[ ) in {x 2 > 0}. 
To conclude the proof, note that by the comparison principle, ( |2.3[ ) and ( |2.2[ ), we have 

C»> " (t ' " C»> * = 2 exp / (1 " dT ) " 2 exp ("^) ' 

if s G (0, 1 /2) and t G [0, 1). Thus, if the solution of (|TTT|)— ([TT2]) obeyed the modulus of continuity 
p for all t G [0, 1), then we would have p(0) = lim^i p(2 — 2t) > 0, a contradiction. □ 

3. The case s g (0, 1) with supercritical drift 

The main ideas for finite time loss of continuity of solutions to ( ] 1 . 1| ) — ( |TT2| ) in the case of bounded 
drifts and s < 1/2 were presented in Section [2j Here we extend those arguments to treat all values 
of s G (0, 1) and corresponding supercritical drifts. For a G (—1, 1) we denote 

{C a (R 2 ) for a G (0,1), 
L°°(M 2 ) fora = 0, (3.1) 
L 2 /H(r2) forae(-l,0). 

In view of the natural scaling of the equations, for any a G (—1, 1 — 2s) the above Banach space 
X a is supercritical. 

Theorem 3.1 (Finite time blow-up in supercritical regime). Let s G (0, 1), a' G (—1, 1 — 2s), 
and let p be any modulus of continuity. Then there exist a smooth divergence-free time-independent 
vector field u with \\u\\ Xa i < 1 and a smooth function 9 with ||$o||c 2 (r 2 ) ^ 1 such that the smooth 
solution of < \ 1 . 1 [ )— ( fTT2| ) breaks the modulus p infinite time, bounded above independently of p. 

Notice that this result contains Theorems | 1 . 1 | and [L2 It does not directly cover the case p G [1,2] 
in Theorem 1.2 but this follows from the case p G (2, 2/(2s — 1)). This is because our vector fields 
will be supported in -84(0, 0) and because the bound < 1 can De replaced by < c 

for any c > due to M in ( |3.4[ ) being arbitrarily large. 

Let us now turn to the construction of the drift u. For a G (—1, 1), define the stream function 

\x\ — x 2 \ 1+a — \xi + x 2 \ 1+a ( X\ 



1pa,s,p( X ^ X 2) = WT\ 1 \ K ~ ^(N) ( 3 - 2 ) 

2(1 + a) \ X 2J 
where the last two factors are smooth cutoff functions (k is also even) designed to remove singu- 
larities at \xi\ = \x 2 \ and |x| = 0. Specifically, k, p a , s , P £ Co°(^) are such that 

X[-l/2,l/2] < « < X[-2/3,2/3] 
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and 

with k even, e a , s , P > to chosen later, and also 



IaO*)I < o 0-3) 



(which is possible for any choice of > 0). Here k removes the singularities at \x\\ = \x 2 \ 



(except the origin) but does not alter the fraction in <\3.2) on the union of two cones given by 
\xi | < \ \%i 1- The subsolution we will use, similar to that in ( |2.2[ ), will be supported in this set, 
so k will not affect the argument. Likewise, p. a , s ,p removes the singularity at the origin and will 
make our drift compactly supported. Since the subsolution will also be supported in the annulus 
B 2 (0) \ B 2£a (0) at all times until p is broken, fJ, a ,s,p will also not affect the argument from the 
previous section. In fact, £ Q , s ,p will be chosen so that the modulus p will be broken before the 
support of the subsolution reaches B 2£a (0). 
We now let, for some M > 0, 

Ua,s,p = J^^ ± ^a,s,p, (3.4) 



which is smooth and compactly supported because so is ip a ,s, P - Then ( |3.3[ ) ensures that for any 
given a G [0,1) there is M > such that ||w aiSiP ||x a < 1 for all s,p, while for any given 
a E (—1, 0) and any a' E (—1, a), there is M > such that s p || Xa ' < 1 for all s, p (the latter 
because if a < 0, then IV"^^^^)! < C^xl" with C a independent of s, p). If now s E (0, 1) 
and a' E (—1,1 — 2s) are given, we either pick a = a' and the M associated with a (if a' E [0, 1)) 
or pick some a E (a', min{l — 2s, 0}) and the M associated with a, a' (if a' E (—1, 0)). In either 



case, the drift from ( |3.4| ) will satisfy H^a^Hx"' — 1 f° r an Y P ( an d M = M ayS is independet of p). 




Figure 1. Velocity field for Figure 2. Velocity field for 

a = -1/2. " a = 1/2. 
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We also note the explicit formula 

u a ,s, P ( x i> x *) = Sg ^ 2 ^ (~\xi - x 2 \ a + \xi + x 2 \ a , -\x l - x 2 \ a - \x t + x 2 \ a ) (3.5) 
for all x G C atStP = {i 6 I 2 | \xi\ < \\x 2 \ and \x\ G [2s aiSiP , 2]}. In particular, we have 

u a , s , p (0,x 2 ) = M~ 1 (0,-x 2 <) (3.6) 

forx 2 G [2e ajS)P ,2]. 

Finally, note that since u a ^ p is smooth and compactly supported, the solution 9 of ( ] 1 . 1 1 ) — ( |TT2] ) 
will be smooth and decaying at spatial infinity uniformly in t G [0, T] for any T > as long as 
0q is smooth and decaying at infinity. We will assume this from now on. In what follows, we will 
also drop the subscripts and denote u = u a , s ,p, £ = £ a,s, P , and C = C a>s>p . 

Lemma 3.2 (Symmetry). If do is odd in x 2 , then so is the solution 9{t,-) of ( | 1 . 1 [ ) — < [TT2] ) for t > 0. 

Proof. This follows from uniqueness of the solution and from u(Rx) = Ru(x), which is due to 
^a,s,p being odd in x 2 . Here and throughout the paper we denote R(xi, x 2 ) = (xi, —x 2 ). □ 

The symmetry of u also shows that if an odd-in-X2 function 9 is a sub(super)-solution of ( |1.1| ) 
in the upper half-plane {x 2 > 0}, then it is a super(sub)-solution in the lower half-plane. In view 
of Lemma [X2| it is therefore natural to compare odd-in-x 2 solutions with odd-in-x 2 functions that 
are sub(super)-solutions on {x 2 > 0} only. From now on we will call such functions odd-in-x 2 
sub(super)-solutions, and for them we have: 

Lemma 3.3 (Comparison principle). If 9 is smooth, decaying at spatial infinity uniformly in 
t G [0, T] for any T > 0, an odd-in-x 2 super- solution of ( |1.1[ ), and satisfies sgn(x 2 )^(0, x) > 0, 
then for t > we have sga(x 2 )9(t, x) > 0. 

Proof. Due to oddness in x 2 and smoothness, we have 9(t, x\, 0) = for t > 0. Assume, towards 
contradiction, that for some T > we have 

inf 9(t, xi, x 2 ) < 0. 

a;ieM,i2>0,te[0,T] 

Since 9 is smooth, odd in x 2 , decays at spatial infinity uniformly in t G [0, T], and ^(0, ■) > 0, the 
above infimum must be attained at some point (t, x), with t G (0, T] and x 2 > 0. We then must 
have d t 9(t, x) < and V0(t, x) = 0, hence (-A) s 9(t, x)>0. 
On the other hand, oddness in x 2 yields 

Jr2 \x - y\ 2+2s 



{V2>0} 



9(t,x)-9{t,y) 9{t,x) + 9{t,yY 
\x-y\ 2 + 2s \x-Ry\ 2 + 2s ' ' U 



(recall that Ry = (y x , —y 2 )). For any y in the upper half-plane we have \x — y\ < \x — Ry\ because 
x 2 > 0, and since the minimum of 9 occurs at x, we also know 9(t, x) — 9(t, y) < 0. Hence 

9(t, x) - 9(t, y) , 9(t, x) + 9(t, y) < 9(t, x) - 9(t, y) + 0(t, x) + 9(t, y) 26 (t, x) 



\x-y\ 2 + 2s \x-Ry\ 2+2s ~ \x - Ry\ 2+2s \x - Ry\ 2+2s \x - Ry\ 2+2s 

and therefore (—A) s 9(t, x) < 2c s 9(t, x) Jr y2>0 \ \x — Ry\~ 2 ~ 2s dy < 0, a contradiction. □ 
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The above lemma will be applied to the difference between an odd-in-x 2 solution 9 and an odd- 
in-x 2 sub-solution 9 as of so that 9 — 9 a s is an odd-in-x 2 super-solution. If the latter is 

non-negative in the upper half-plane at t = 0, it follows from Lemma [33] that at all later times, the 
oscillation of 9 over any disc centered at the origin is no less than the oscillation of 9 a _ s over the 
same disc. 

We now turn to the construction of 9 a , s , in the spirit of the argument in Section]!] Let 77 G C°° (R) 
be an even function, supported on [—1, 1], positive on (—1, 1), non-increasing on IR + , and with 
77(0) = 1. For some r a G (0, 1/4) (to be chosen later), we let 

<j> a (x 1 ,X 2 ) = 77(1(^1,^2 - l)ka X ) - V(\( X U X 2 + I)!?"" 1 ) (3.7) 

be a smooth approximation of XB ra (o,i) — XB ra {o,-i)- As in Section[2j we shall use <\> a to build 9 a>s , 
but first we give a measure of the effect of (— A) s on <\> a . 



Lemma 3.4 (Control of dissipation). Let a G (—1, 1), r a G (0, 1/4), and (f) a from ( ]3.7[ ). Then 
there exists c atS > such that 

(-A) s Q (x) < c a ^ a (x) (3.8) 
holds in the upper half-plane {x 2 > 0}. The constant c Q>s depends only on r a , s, and r\. 



Proof of Lemma 3.4 We start with x G {x 2 > 0} \ B ra (Q : 1). Due to the definition of 77, we have 
<j>a( x ) — 0' an d using oddness in x 2 we obtain 

(-A)v«W = -c,/ t5 ^L-d» 

Jr 2 F — 1/1 

1^(0,1) |x-7/| 2 + 2s JB ra {o-i) \x-y\ 2+2s 

„, M ' (I( » - 1)|r "» - |^^) (3.9) 

Now we notice that for x,y in the upper half-plane we have \x — y\ < \x — Ry\, and hence the 
integrand in ( [3T9] ) is positive. It follows that (-A) s (fi a (x) < for all x G {x 2 > 0} \ £v„(0, 1). 
Since dB ra (0, 1) is compact and a G C°°, it follows that there is c\ = c\(r a , s, 77) > such that 

(-A) s Q (x) < -ci < 

for all x G dB ra (0, 1). Hence there exists 5 = 5(r a , s, 77) G (0, r a /2) such that (— A) s (J) a (x) < 
for all x G {x 2 > 0} \ B ra ^s(0, 1), implying ( |3.8[ ) for these x (with any c a)S ). 

It is left to verify ( |3.8[ ) for all ,x G S ra _5(0, 1). in view of the monotonicity of 77, for each such 
x we have that (f> a (x) > 77(1 — ^r^ 1 ) > 0. Therefore, for these x we have 

(-A) s Q (x) < c 2 r~ 2s 1 1 77 1 1 c * < c QiS 7?(l - 5r~ l ) < c a>s (j) a {x) 

for a suitably chosen c a;S > (depending on r a , s, S, 77), which concludes the proof. □ 

From homogeneity of the kernel associated with the fractional Laplacian (— A) s and from 
Lemma [34 we directly have the following. 

Corollary 3.5 (Effect of dissipation under rescaling). For A > let <f> a> \(x) = <f> a (x/X). Then 

(-A) s a , A (x) < c QiS A- 2s a , A (x) (3.10) 

holds on {x 2 > 0}. 
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Having established Lemma 3.4 and Corollary 3.5, it is clear from the argument in Section [2] 
how we will proceed. We will build an odd-in-x 2 sub- solution to ( |1.1[ ) by adding an appropriate 
exponentially decaying prefactor (chosen using p.!0[ )) to an odd-in-x 2 sub-solution of the pure 
transport equation (d t + u • V)0(t, x) = with initial condition <p a . This sub-solution will be a 
time-dependent rescaling of Q , that will depend on the vector field u = u atS:P defined in < |3.4| ). 
Note that it is easier to build a sub-solution then to find the actual solution because the disks on 
which <p a is supported deform when transported by u. 

Let T = T a>StP = M(l - (Ae) l - a )/(l - a) (recall that e = e a ^ p ) and for [0, T] let 

Za (t) = (1 - (1 - a) AT 1 *) 1/(1-00 . (3.11) 



Recalling ( |3.6[ ), we find that z a is the solution of the initial value problem 

z{t) = u 2 (0, z{t)) = -M~ l z{t) a , z(0) = 1, (3.12) 
because it is decreasing and z a {T) = As > 2e. 

Lemma 3.6 (Subsolution of the transport equation). The function <j) a (t,x) = <fi a (x/ 1 z a (t)) is 

an odd-in-x 2 sub-solution of (d t + u a ^ S)P ■ V)0 = on the time interval (0,T a s p ), provided 
r a G (0, 1/4) is sufficiently small. 



Proof of Lemma 3.6 By the definition of 4> a it suffices to prove 

d t <f) + u ■ V<p < 

where 

4>(t,Xi,X 2 ) = V 77^ 

V r a z a (t) 

is the component of (j) a (x / z a (t)) which is supported in the upper half-plane. Let 

v(t,Xi,x 2 ) = —M~ 1 z a (t) a ~ 1 (x 1 ,x 2 ) 

be the vector field that advects 0, so that we have d t (p + v ■ V0 = 0. Hence, in order to prove the 
lemma it suffices to show 

(v(t, xi,x 2 ) - u(xi, x 2 )) ■ V0(t, Xi, x 2 ) > 
which, since r) is non-increasing on IR + , is equivalent to 

(u{xi,x 2 ) + M- 1 2 Q (t) a - 1 (xi,x 2 )) •(x 1) x 2 -z a {t)) >0. (3.13) 

It suffices to prove ( |3. 13| ) for t G [0, T] and x\ + (x 2 — z a (t)) 2 < (r a z a (t)) 2 (notice that such 
(xx,x 2 ) then lie in the domain C because r a < 1/4 and z(t) G [4e, 1]). By a-homogeneity of u in 
C (see p.5[ )) and after scaling by z^t)' 1 , it is sufficient to verify that 

(Mu(x 1 ,x 2 ) + (xi,x 2 )) ■ (x h x 2 - 1) > (3.14) 

holds for x\ + (x 2 — l) 2 < r 2 a . Since p.5[ ) holds for these x, the left-hand side of ( |3.14| ) there 
equals 

f(x%, x 2 ) = --(x 2 — xi) a (x 2 - 1 + xi) - -(x 2 + x 1 ) a (x 2 -l — xi)+ x 2 {x 2 - 1) + x\. 

After an elementary computation we obtain 

/(0, 1) = / Xl (0, 1) = f X2 (0, 1) = f XlX2 (0, 1) = 0, (3.15) 
f XlXl (0, 1) = 2(1 + a) > 0, /^(0, 1) = 2(1 - a) > 0, (3.16) 
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and therefore (0, 1) is a local minimum of / for any a G (—1, 1). Since / is C 2 near (0, 1), there 
exists r a G (0, 1/4) such that f(xi,x 2 ) > for x\ + (x 2 — l) 2 < r 2 a , which proves p.!4[ ), and 
hence the lemma. □ 



Proof of Theorem \3J\ We now let 



c a>s M{\ -a- 2s)- 1 (1 - z a (t) l - a - 2s ) a + 2s < 1, 
—c atS M In z a (t) a + 2s = l, 

(3.17) 



where c ayS is from ( |3.8[ ), and t G [0, T]. Finally, we define a rescaled modulated version of <j) a : 

0a,s(t,x) = Cexp (-7a, s (t)) 0a ' (3.18) 

where c = ||</> a ||^. Again, ajS is a function supported at any time t on two discs whose centers 
(0, ±z a (t)) are transported towards the origin with the drift u and whose radii are r a z a (t), that is, 
they shrink at the same rate as they approach the origin. 

Lemma 3.6 and Corollary 3.5 together show that 9 a , s is an odd-in-x 2 sub-solution of ( |1.1| ) for 
t G [0, T] because then 

(d t + u ■ V + (-A) s ) 9 a , s < %, s (t)9 a , s + c r , s z a (t)- 2s 9 atS = 

holds on {x2 > 0}. If 9 is the odd-in-X2 solution with initial condition 9 = c<f) a (= 9 a)S (0, •)), 
Lemma 3.3 shows that sgn(x 2 )6 l > sgn(x 2 )9 a:S for t G [0,T]. We now only need to pick e = 
£ a,s, P > so that 9 a s (and thus also 9) breaks p at the origin at time T = T a s p . This will be 
possible because if we picked e = and thus T = M/(l — a) (recall that a + 2s < 1, so 
z a (M/(l — a)) = 0), then we would obtain 7 QjS (T) < oo, so 9 would become discontinuous in 
finite time. 

Of course, we will need e > to ensure smoothness of the drift u and the solution 9. We let 

£ = £a,s, P > be such that 

p{8e) < 2cexp(-7« iS (M/(l - a))) = 2cexp (-c a , a (l - a - 2s)" 1 ) (> 0) (3.19) 
and consider the corresponding T = p < M/ (1 — a). Then we immediately obtain 

osc 0(T,-)> osc # QS (T, •) > 2cexp(- 7as (T)) > 2cexp (-c Q Jl - a - 2s)" 1 ) > p(8e) 

because 9, 9 aiS are smooth. Thus p is broken at time T < M/(l — a), with M independet of p. □ 

Remark 3.7 (On bounds for the critical case). In the critical case a + 2s = 1, the above proof 
fails because if T is such that z a (T) = 0, then 7 QjS (T) = oo. However, ( |3.17| ) does yield 



osc 9(t, ■) > 2cexp (— 7q!)S (£)) = 2c2 a (t) 



C„..M 



for t G [0,Tq, s p ]. Thus, the solution 9 cannot have Holder modulus of continuity C 13 for any 
(3 > c ajS M, so it becomes less regular as M decreases. In particular, for s = 1/2 and a = 0, we 
cannot have a Holder modulus better than C||w||^i, with a universal C > 0. 

The procedure described in this section may be summarized in one abstract theorem, which we 
shall use later in Section [4] to obtain blow-up in the case s = 1, and in Section [7] for the case of a 
nonlocal operator £ which generalizes the fractional Laplacian, and is slightly supercritical. 
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Theorem 3.8 (An abstract loss of regularity result). Assume that the drift-diffusion equation 

e t + u ■ w + ce = o (3.20) 

on (0, T) x IR 2 , with u a continuous divergence-free vector field with u(t,Rx) = Ru(t,x) and 
C a dissipative linear operator acting on x, satisfies Lemmas \3.2\ and \3.3\ Further assume that 



Ui(t,0,x 2 ) = and sgn(x 2 )u 2 (t 1 0, x 2 ) < for x 2 ^ 0. Let r G (0,1), rj be a smooth bump 
function as above, define 

4>(xx,x 2 ) = r)(\(x 1 ,x 2 - 1) |r _1 ) - r)(\(xi,x 2 + 1) |r _1 ), 

and let z solve the ODE 

z(t)=u 2 (t,0,z(t)), z(0) = 1 (3.21) 
on (0, T). Assume that for A > and (fi\(x) = <p(x/X) there exists H(X) > such that 

C<p x < H{\)4> x (3.22) 
on R x ]R + . Finally, assume that for each t 6 (0, T), 

xi,x 2 ) - u 2 (t,0, z(t))z(ty 1 (xi, x 2 )) ■ (xx,x 2 - z(t)) > (3.23) 

holds on the disc x\ + (x 2 — z{t)) 2 < rz(t). If 9 is the solution of ( |3.20[ ) with initial condition 
8(0, •) = c<f)for some c > 0, then 9 breaks at time T any modulus of continuity p with 



p{2z(T)) < 2cexp - / H(z(t)) dt 



(3.24) 



Notice that if z(T) = (which can only happen if u is not uniformly Lipschitz in x) and 
f H(z(t)) dt < oo, then 9 breaks any modulus of continuity at time T, so we have finite time 
blow-up inC(M 2 ). 

It is clear that Theorem 3. 1 is a particular case of Theorem 3.8 with C = (— A) s and u = m QjSiP 
from p.4[ ). Lemma 3.4 and Corollary 3.5 merely verify ( |3.22| ), while ( |3.23[ ) is verified in the proof 
of Lemma 3.6 Finally, super-criticality a + 2s < 1 guarantees condition ( |3.24| ) for any given 
modulus p and an appropriate T = T a:S:P . 



Proof of Theorem \3H\ This follows as above because 4>(t,x) 



(x/z(t)) is again an odd-in-x 2 



sub-solution of (d t + u ■ V)0 = for t £ (0, T) (see the proof of Lemma 3.6 up to ( [3T3 >), and so 



9(t, x) = cexp 



H{z{r))dT 



x 



ccj). Hence again 



(3.25) 



is an odd-in-a;2 sub-solution of ( |3.20[ ) for t G (0, T) with ^(0, • 

osc 9(T,-)> osc 9(T, ■) > 2cexp (— [ H(z(r)) dr) > p(2z(T)\ 

B Z(T )(0) B Z(T) (0) V JO J 

so 9 breaks p at time T. 



□ 
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4. Finite time blow-up for non-autonomous drift-diffusion 



In this section we prove Theorem 1.3 The argument in Section [3 fails for s — 1 when the drift 
velocity u is time-independent. In fact, as we shall prove in Section 6] that the result is simply not 
true in this case: the solution has a very weak modulus of continuity, globally in time, even if the 
time-independent drift is supercritical. In this section we prove that by allowing the drift velocity 
to depend on time, finite time loss of regularity can still be obtained in the supercritical regime. 
We recall that for s = 1 the critical space is L 2 (IR 2 ). 



Proof of Theorem \L3\ The proof is an application of Theorem 3.8 We first let q E (p, 2) and let 
z(t) be the solution of 



z(t) = -M- 1 z(t)- 2/q 
with M > to be chosen later. That is, 

q + 2 



z(0) 



z(t) 



1 



M~H 



Q 



(4.1) 



(4.2) 



for all t G [0, t q ], where t q = Mq/(q + 2) is the time when z reaches 0. 
Next, we consider the vector field 



u(x) 



Fi - x 2\ - Fl + X 2\ 



%1 
•i'2 



MM) 



(4.3) 



where k, \i are smooth with k as in ( |3.2| ) and 

X[l/2,2] < [l < X[l/3,3]- 

That is, u is a smooth and compactly supported version of the vector field in Section [2] 
Finally, we let u be any smooth time-dependent divergence-free vector field with 



u(t, x) 



Mz(t) 2 /i \z(t) 



(4.4) 



for t G [0, T p ] and reasonably behaved for t > T p , where T p e (0, t q ) will be chosen later. We now 
pick M > large enough so that || sup t \u\\\lp < 1. This is possible because 



sup \u(t, x)\ < 

t<T a 



C 



M\x\ 2 /i 

(for some constant C), and the left hand side vanishes for \x\ > 3. 



Then u and C = —A satisfy all the hypotheses of Theorem 3.8 with T = T p , any r G (0, 1/4), 
as well as for instance 77 ( | x |) = exp(l — (l — |a;|) _1 )x[o,i](l :r l) an ^ c > such that then Hc^Hc-z < 1. 
In particular, comparison and oddness in x 2 follow; p.21[ ) holds because u(0, 1) = —(0, 1) and 
d3^3) is just ((23) (both fo rt < T p ). 

It remains to check that ( |3.22[ ) holds with H(X) = C\~ 2 for some C > 0. By scaling, we only 
need to check it for A = 1. Since we are in the case s = 1, which is local, only one of the terms in 
the definition of <p affects A0 in the upper half plane: 

A(j)(x) = Ar}(r~ 1 \(x 1 ,x 2 - 1)|). 
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We start checking ( |3.22| ) for x sufficiently close to dB r (0, 1). For x with d = r 1 \(xi,x 2 — 1] 
close to 1, we have 



-A<f>(x) 



V"(d) 
1 



V'(d) 



[l-dy (l-df d( 



T^m) exp(1 ~ (1 ~ dyl) - °- 



The last inequality holds for d > 1 — S if 5 is a small constant, since the first term (1 — d) 4 
controls the other two. For x such that r~ l \ (xi, x 2 — 1) | < 1 — 5 , we have 

-A(j)(x) < r- 2 ||Ar7|| L oo < Cr}(l - 8 ) < Ccj)(x) 

by the monotonicity of r\, if C is chosen sufficiently large. 

Since H(z(t)) = Cz(t)~ 2 , the right-hand side of ( |3.24| ) is strictly positive even if T is replaced 
by t g (because —2q/(q + 2) > —1). Thus one only needs to choose T p sufficiently close to t q so 
that ( |3.24[ ) holds. The proof is finished. □ 



5. Uniqueness of particle trajectories for divergence-free time-independent 

drifts in two dimensions 

In this section we show that when the velocity of a transport equation is both divergence-free 
and time-independent, a seemingly unexpected regularity result for the associated Lagrangian tra- 
jectories holds. 

Given a time-independent vector field u which is divergence-free and does not vanish on a two- 
dimensional domain, we let X(a, t) be the flow map associated to u, i.e. the solution of the ODE 

dX ^^ =u{X{a,t)), X{a,0)=a. (5.1) 

Although u is not Lipschitz, since u is divergence-free, we will show that the flow map X(a, t) is 
unique, as long as \u\ stays away from 0. More precisely, we have the following: 

Theorem 5.1 (Uniqueness of particle trajectories). Let u : IR 2 — > R 2 be a divergence-free con- 
tinuous function such that \u\ ^ on a domain T> C M 2 . For any a G T>, the solution of (5J_) is 
unique (forward and backward in time ) as long as it stays in V. 



Proof of Theorem \5J\ Existence of solutions is standard because u is continuous (Peano's theo- 
rem). We prove uniqueness by contradiction. Modulo changing the direction of time (i.e., replac- 
ing u by —u), we can assume that there exist a ^ b E V and a minimal time T > such that 
X(a, T) = X(b, T). Here T is assumed to be small enough so that X(a, t),X(b, t) 6 V for all 
t E [0, T]. Let u = u(X(a,T)) = u(X(b,T)). Since u is continuous, for small t > we have 
X(a, T — t) = X(a, T) — u t + o(t). Therefore, for any small e > 0, we can find £i and t 2 small 
such that 

(X(a, T) - X(a, T - h)) ■ u = (X(b, T) - X(b, T - t 2 )) ■ u = e. 



We define the paths 
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71 = {XfaT-h + s): s G [0,^]} 

72 = {X(6,T- S ): se [0,t 2 ]} 73 

73 = {sX(b, t 2 ) + (1 - s)X(a, h):se [0, 1]}. 



7i 



n 



X(a,T) 



72 



X(b,T-t 2 



Due to the minimality of T, it follows that 7 = 71 U 72 U 73 is a piecewise smooth simple closed 
curve (let us orient it positively). Call the region enclosed by 7 and let n be the outward unit 
normal vector to dVl. Note that on 71 U 72 we have u ■ n = 0, and hence, applying Green's theorem 
in Vt yields 







V • u 



In the last inequality we use that on 73, n 
then u is close to u along 73. 



u ■ n = I u ■ n < 0. 

•J 73 

Uq/ I uq I and, since e is small and u is continuous 

□ 



Remark 5.2. Note that Theorem 5.1 immediately implies that the transport equation 

d t 9 + u ■ V0 = 0, 

has a continuous solution 9 if the initial data 9 is continuous and the drift u is time-independent, 
continuous, divergence-free, and non-vanishing. 



6. A MODULUS OF CONTINUITY FOR AUTONOMOUS DRIFT-DIFFUSION 



In this section we prove Theorem 1.4 That is, we consider the equation 

d t 9 + u(x) ■V6-A6 = 
0(0, •) = #0 



(6.1) 
(6.2) 



in the full space 



x 



2 , with a divergence-free vector field u E L\ oc . It is important here that u 



is time-independent, and it will be used in the proof that d t 9 satisfies the same equation. 

We say that such u is divergence-free (in the distributional sense) if L 2 u(x) ■ V(j)(x)dx = for 



each 



eC °°( 



), and 6 eL c 



is a distributional solution of ( ]6. 1 [ ) — ( f6T2| ) if 



9(t, x)((f) t (t, x) + u(x) ■ V0(t, x) + A(p(t, x))dtdx 



9 (x)(f)(0,x)dx (6.3) 



for each <fi G C^°(IR 3 ). Note that it is enough to prove the a priori estimate in Theorem 1.4 for 
smooth u and smooth classical solutions. Indeed, the result of the theorem then follows by a 
standard approximation of the drift u in the L\ oc norm by smooth divergence-free drifts u n , for 
instance, via a mollification. The uniform a priori modulus of continuity for the associated smooth 



solutions 9 n (spatial from ( 1.3 ) and temporal from the last bound in Lemma |6T2[) and the maximum 



principle for 9 n show that 9 n converge (along a subsequence) locally uniformly in 



to some 



9 e L°°(R + x 
t — > 0, since all the (9 n ,u r 
it follows that the right-hand sides of ( |6.3[ ) for (9 n , u n ) converge to the right-hand side of ( |6.3| ) 



satisfying the same modulus of continuity. Despite this modulus blowing up as 

is compactly supported, 



are uniformly bounded in L°° x L\ oc and 
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for (9,u). Thus 9 is a distributional solution of (|6.1|)-(|6.2|). Note that we do not claim that the 



weak solution of the problem with u G L loc but not smooth is unique, and there might be other 
distributional solutions, which could be discontinuous if u is not in L 2 . If u is not at least in L\ oc , 
the definition of divergence-free fields and distributional solutions is questionable. 

In the rest of this section we will assume that u is smooth and hence so is 9. The idea of the proof 
is as follows. It was noted in OTTl that any function in IR 2 which satisfies the maximum principle 
on every disk and belongs to H 1 has to be continuous, with a logarithmic modulus of continuity 
depending on its H 1 norm only. The H 1 estimate, which comes from the energy inequality, is 
a critical quantity in terms of continuity of the solution, since any estimate in H 1+£ (K 2 ) would 
imply a Holder modulus of continuity from classical Sobolev embeddings. It is well known that 
a function in i/ 1 (IR 2 ) is not necessarily continuous. However, if the function solves some elliptic 
PDE, the maximum principle can be used to bridge that gap and obtain a logarithmic modulus of 
continuity. 

For the case of parabolic equations, we have two difficulties in carrying out this scheme. The first 
one is that the energy inequality only gives us an estimate in L 2 (if 1 ), which is far from implying 
continuity. However, by differentiating the equation in time (and using that u depends only on 
x) we can get that the solution is actually bounded in C(H l ) and Lipschitz in time. The other 
difficulty is more severe, and it is that the parabolic maximum principle is quite different from the 
elliptic one. For example, if we fix any function / G H l (Jst 2 ) and extend it as constant in time, 
this function belongs to C(H l ), is certainly Lipschitz in time, and satisfies the parabolic maximum 
principle without necessarily being continuous. In order to overcome this second difficulty we 
prove that at each fixed time t, the solution approximately satisfies the elliptic maximum principle 



(Lemma 6.8 ). In order to obtain this approximate maximum principle at each t, we need to use the 



equation again, applying a non-trivial result due to J. Nash (see Theorem 6.1 and Corollary |6.6 ) 



6. 1 . Energy estimates. We start by discussing some estimates on the fundamental solution G of 



( |6.1[ ). Recall that G is by definition the solution of the equation 

d t G(t, x, y) + u(x) ■ V x G(t, x, y) - A x G(t, x, y) = 0, 
with lim^o G{t,-,y) = 5 y . With respect to y, it solves the dual equation 

d t G{t, x, y) - u{y) ■ V y G{t, x, y) - A y G(t, x, y) = 0, 
with lim t _n) G(t, x, •) = 5 X . 



The fundamental solution is used to compute the solutions of ( |6.1| ) via the formula 

9{t,x) = [ G{t,x,y)9 {y)dy. 



It is a simple consequence of the maximum principle that G is nonnegative and J G(t, x, y) dy = 1 
for all x and t. The following is a result by J. Nash on the size of the fundamental solution. The 
remarkable fact of this estimate is that it is independent of u, as long as it is divergence-free. 

Theorem 6.1 (J. Nash QUI ). Ifu is divergence-free in IR 2 , then the fundamental solution G(t, x, y) 



of ( |6.1| ) satisfies the pointwise bound 

< G{t,x,y) < C/t, 

where C is a universal constant (independent ofu). 

The theorem was not stated explicitly with these assumptions by Nash, therefore we include the 
proof below. It is a straight-forward adaptation of the proof in ||30l to our setting. 
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Proof. We compute the evolution of the L 2 norm of G. For any fixed y, let E(t) = \\G(t, •, y) ||| 2 
(which is finite for t > 0). We have that 

E'(t) = 2 [ G(t,x,y)(-u- V X G + A x G)dx = -2 / iV^G^x^pdx < -C^Efr) 2 . 
Jr 2 Jr 2 

for some universal C > 0. In the last inequality we used the interpolation inequality (Nash's 
inequality) 



T d < 11/11* \\f\\t d , 



together with the fact that \\G(t, -, y)||z,i = 1 for all t and y. 

The ODE £"(t) < -E(t) 2 /C implies that £7(t) < C/t. Thus, we have shown a universal bound 
•, y)\\ 2 L 2 < C/t. Since the estimate does not depend on u (in particular it is the same if we 
replace u by — u) we also have \\G(t, x, -)\\ 2 L 2 < C/t. 

The pointwise estimate now follows from the energy estimate via the semigroup property of G: 



G(t, x,y)= / G(t/2, x } z)G(t/2, z, y) dz < \\G(t/2, x, •) \\ L 2 \\G(t/2, •, y) \\ L , < C/t. 
Jr 2 

The proof is finished. □ 



Lemma 6.2. Let 9 be a smooth solution of (6.1 ) in R + x IR . Then 



\\Q{t, OIUh* 2 ) - II^oIUhr 2 ) 

II#(V)||l°°(R 2 ) — ||#o||l°°(R 2 ) 

-)||il(K2) < (l + |k||il oc (R2)J ||6 I o||h/4,1(R2) 

l^t^(^) -)IU°°(M 2 ) < ||#o||c 2 (R 2 ) + Ct IMU/ (R 2 ) || ^0 || VF 4 ' 1 



for some universal C > 0. 



Proof. The first two estimates are just the maximum principle for parabolic equations, which also 
holds for all LP norms with p < oo because u is divergence-free. 
Since the drift is time independent, d t 9 solves the same equation 

d t (d t 6) + u-V(d t 6)-A(d t 6) = 0. 

The third estimate follows now directly from d t 9(0, •) = — u ■ V6* + A# and the maximum 
principle, using also 

\\u -V^lln < \\u\\li c ^2 

nez 2 

<c\\u\\ Ll (||V0 o |Ui + ||V 4 0o|Ui) <c\\u\\ri \\e \\ w *,i. (6.4) 

loc v / loc 
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Note that if instead of u E L\ oc we assume that u E L 1 , then one may lower the regularity 
assumption on the initial data to 9 E C 2 D W 2 ' 1 . Finally, if G is the fundamental solution, then 



d t 9(t,x) = I G(t,x,y)d t 9(0,y)dy 

G(t, x, y)(-u(y) ■ V9 (y) + A9 (y)) dy 



< ||A0„||l<» + j G(t,x,y) \u(y) ■ V9 (y)\ dy 

< \\e \\c* + \\G(t,-,-)\\ L ao\\u-V6 \\iA 

< ||0o||c 2 + C^IMIl, 1 ll^ollw 4 -^ 



loc 



where in the last inequality we applied Theorem 6. 1 , and in the one before we used (6.4 ). □ 



By Lemma 6.2 both 9 and d t 9 are bounded in L 1 fl L°° for positive time, and in particular, they 
are bounded in L 2 for every t > 0. The following lemma can then be applied to obtain estimates 
in space-time norms. 

Lemma 6.3. Let 9 be a smooth solution of ( |6.1| ) in R+ x M 2 . Then 

P\\l?({t Q ,oo),fr) < l|0(*O, OIU 2 
ll3t0|lL i »((to,oo),H 1 ) ^ ll^(*0, •) IU 2 
ll e IL-(( i0l oo),Hi) < C{\\9{t Q ,-)\\ L 2 + ||^(t ,-)IU 2 ) 

for some universal C > 0. 

Proof. The first two estimates are the classical energy estimates for the equations for 9 and d t 9 
respectively. Note that since u is divergence-free, the drift term has no effect on the energy esti- 
mate. The first and second inequalities can be interpreted as that 9 E -ff 1 ((0, oo), H 1 ). The third 
inequality follows therefore from the one dimensional Sobolev embedding H 1 <^-> C 1 ! 2 . □ 

Remark 6.4. The bounds on u in L\ oc and 9 in C 2 fl W 4 ' 1 (alternatively, on« E L 1 and 9q E 
C 2 n W 2 ' 1 ) are only used to obtain estimates on d t 9 in L°°((t , oo); L°°(R 2 )) and on V9 in 
L oo ((t , c>o); L 2 (IR 2 )) for any t > 0. These estimates can also be obtained in terms of, for in- 
stance, 

||6> ||l 2 = Ai < oo and ||it • V^o — A9 \\ L 2 = A 2 < oo. 

Indeed, from d t 9(t,x) = J G(t, x,y)d t 9(0, y) dy and \\G(t, x, ■) || £a < Cr 1 ' 2 , it follows that 
\\d t 9(t, ^Wl^ < Ct _1 / 2 v4 2 . Also, the maximum principle for the equation obeyed by d t 9 implies 
\\dt9(t, ■) ||l 2 < A 2 , and from the energy inequality we get an estimate on d t 9 in L 2 ((t , oo); H 1 ^ 2 )). 



The bounds in Lemmas 6.2 and 6.3 follow, and they depend on the (smooth) drift u only implicitly 
through A 2 . 

6.2. A maximum principle on time slices. In this sub-section we denote B r = B r (0) C M 2 . 

Lemma 6.5. T/zere exists h > swc/j that for any smooth divergence-free u on B\, the solution of 

d t 4> + u ■ V0 — A0 = m (-/i, 0) x Bi 
4>(—h, •) = 1 in -Bi 

= in (-/i, 0) x dB 1 

satisfies sup Bl 0(0, ■) < 1/2. 
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Proof. The function <p is nonnegative, so if we extend it by zero outside of Bi for all times t (and 
extend u in any way we like), we obtain a sub-solution of the equation ( |6.1[ ) in space (—h, 0) x M 2 . 
Therefore, by the comparison principle, it has to stay below the actual solution with the same initial 
values at {— h} x R 2 : 

4(P, x)< [G(h,x,y)x Bl (y)dy<ch-\ 



where we applied Theorem 6T]in the last inequality. So h = 2C works. □ 



6.5 



□ 



Corollary 6.6. With hfrom Lemma 6.5 we have that for any a > b, the solution of 

d t <p + u ■ V<f) - A(p = in {-hr 2 , 0) x B r 
<f)(—hr 2 , •) < a in B r 

<f><b in(-hr 2 ,0) x dB r 

satisfies sup Br 0(0, ■) < (a + b)/2. 

Proof. Compare <fr with b + (a — b)cj)(r 2 t, rx)/2, where cf> is the function from Lemma 

We recall the maximum principle for parabolic equations: 
Proposition 6.7. For any parabolic cylinder Q = (—h,0) x B r , we have that 

sup 9 = sup 9, 

Q d p Q 

where d p Q denotes the parabolic boundary 

d p Q = ({-h} x B r ) U ((-h,0) x dB r ) . 

The following lemma gives a relation on each time slice which approximates the maximum 
principle for elliptic equations. 

Lemma 6.8 (A maximum principle on time slices). Let h be the constant from Lemma |63] and 
assume that 9 is a solution of ( |6.1[ ) in [— T, 0] x Bn^, with d t 9 bounded. Then for any R < \/T/h 

we have 

sup 0(0,-) < sup 0(0,-) + 2\\d t 9\\ L oohR 2 , 

Br dB R 

mf 0(0, •) > inf 9(0, ■) - 2\\d t 9\\ L o*hR 2 . 

Br 8Br 

Proof. We only prove the first claim, the second follows by considering —9 instead of 9. Let a = 
sup Bfl 9(—hR 2 , ■) and b = sup 9BjJ 9(0, ■) + \\9 t \\ L oc,hR 2 . Then obviously b > max(_ fcfi 2^ xaBji 9. 
If a < b, the maximum principle immediately yields sup ™ 9(0, ■) < b so let us assume a > b. 
On one hand, from Corollary 6.6 we have that 

a/a \ ^ a + b / a + sn PdB R 0(P, ■) + \\0th°°hR 2 
sup 9(0, ■) < < 



b r 2 2 

On the other hand, from boundedness of d t 9 we also have 

a < sup 9(0, ■) + \\9 t \\ L00 hR 2 . 

Br 

Combining the last two inequalities, we obtain the first claim. □ 
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6.3. A local modulus of continuity and the proof of Theorem |1.4[ From Lemmas |672| and [63 



we have that the hypotheses of Theorem 1.4 imply that ||#||z,°°, \\d t 9(t, -)IU°°' an d ||#(t, Ollff 1 WQ 
bounded on any time-interval [t , oo) with t > 0. The proof will be finished by using the following 



local result. It applies Lemma 6.8 on each time slice, so that we can now more or less follow the 
idea from the proof in 113111 . 

Theorem 6.9 (Local continuity with supercritical drift). Let u be a smooth time-independent 
divergence-free vector field on -Bi(O) C IR 2 and assume that the solution 9 of ( |6.1[ ) on (t , ti) x 
5i(0) satisfies \\d t 6(t, •) \\ L oo( Bl (o)) < D, || V6»(t, ■) || L 2 (Bl(0)) < E, and ||0(*, 0IU<»(Bi(o)) < ^ 
uniformly in t G (to, 77?en 0(t, •) restricted to .81/2(0) satisfies a modulus of continuity given 
by 



Pt{r) 



C(D + E + F^\og_(t-t ) 



(6.5) 



\J- log r 

/or a// r G (0, 1/2) and t G (to, ti), and some universal C > (wzY/j log_ s = max{— log s, 0}). 



Proof of Theorem [679] Fix x G £>i/ 2 (0) and t > to, and let T = min{t — to, h/4}, with h from 
Lemma 6.5[ Let also c = We now apply Lemma 6.8 to 9 shifted by (t, x) and get for any 



s < ^T/h < 1/2, 



In particular 

We first prove that 



sup 9(t, 

8B s (x) 

inf Bit, 

dB 3 (x) 



osc 9 it, 

8B a (x) 



> sup 9{t, ■) — cs 2 , 

B s (x) 

< inf 9(t, ■) + cs 2 . 

B s (x) 



> osc 9(t, 

B a (x) 



2cs 2 . 



(osc 9 

B r 



osc 9 it, 

B r (x) 



< max < 4cr 



47TE 



V-lOg': 



(6.6) 



(6.7) 



holds for any r < Tjh{< 1/4), after which we proceed to all r < 1/2. To prove ( |6.7[ ) we just 
need to consider the case when oscs r 9 > Acr. Let R = y/r < >jT/h and estimate 

R R 

E 2 > J |W| 2 dx = J J |V0| 2 dads> J J \9 a \ 2 dads, 

B R \B r r 8B S r 8B S 

since | V9\ 2 = 9 2 + 9 2 where 9 a is the tangential derivative and 9 U is the normal one. We rewrite 
the integral on the right using polar coordinates (sd = a), and use Cauchy-Schwartz to obtain 



R 



E > 



h(sa)\ 2 dads > 



- (^osc I 

1T 2 S 8B S 



ds. 



8B X 



Applying estimate (|6.6[), and noticing that 2cs 2 < 2cR 2 = 2cr < \ osc Br 9, we get 



R 



7T 2 E 2 > 



-(osc 9 — 2cs 2 ) 2 ds > (logi? — logr) ( - osc 9 

S B r \ Z B r 



-logr 



22 



LUIS SILVESTRE, VLAD VICOL, AND ANDREJ ZLATOS 



Thus < |6.7| ) holds for r E (0, T/h). On the other hand, we have that 

U°°(B r ) < 2F 



osc < 2 



(6.8) 



holds for r 6 pT//i, 1/2). Thus, in order to prove ( |6.5| ), we only need to combine ( |6.7| ) and ( |6.8| ) 
with the fact that r < (— log r) -1 / 2 for r < 1/2. □ 



The proof of Theorem 1.4 now follows by fixing any x G IR 2 (without loss take x = 0) and 
t > 0, then letting t = t/2 (with ft, from Lemma 6.5 ) and tx = oo. Lemmas 6.2 and [63] imply that 



d < c(i + r 1 \\u\\ Ll j\\e \\ O 2nw^ 

E < C(l + \\u\\ ri + r 1 ||w|| z i )||6 l o|| ( 72niy 4 . 1 

loc loc 



which combined with (6.5 ) implies the estimate ( 1.3 ) of Theorem 1.4 Lastly, by Lemma 6.2 9 is 



x 



uniformly Lipschitz in time on [to, oo) x IR 2 , for any to > 0, and hence continuous on 

7. LOSS OF REGULARITY IN THE SLIGHTLY SUPERCRITICAL CASE 

Motivated by [fT9l , in this section we address the regularity of solutions to the drift-diffusion 
equation 

d t e + u-ve + ce = o (7.i) 

with smooth initial condition 9(0, x) = 9q(x) and a bounded divergence-free drift u. Here C is 
a nonlocal dissipative operator which is slightly less smoothing than (—A) 1 / 2 . More precisely, 
let m: IR 2 \ {0} — > (0, oo) be a smooth radially symmetric, radially decreasing function, that is 
singular at the origin, decays at infinity, and satisfies the below properties: 



m(r) 



dr < +oo 



1 + r 

is non-decreasing for r e (0, 1). 
We abuse notation and write m(y) = m(\y\) for y EM 2 \ {0}. Condition ( |7.3| ) can be relaxed, to 



'o 

rm(r) 



(7.2) 
(7.3) 



r P m i r 



nonlocal operator 



is non-decreasing on (0, 1) for some (3 < 2. Associated to this function m we define the 



C9(x) = P.V. 



(9(x)-9(x + y)) r ^dy 

\y\ 



(7.4) 



for all smooth functions 9. Condition ( |7.2[ ) is the only essential one, and shows that we may take 
m(r) ~ r~ s for any s E (0,1/2), but also m(r) w r~ 1 (log(2 + l/r))^ 13 for any (3 > 1. We 
informally say that C is less smoothing than (—A) 1 / 2 by at least a logarithm, or that C is slightly 
supercritical with respect to the scaling induced by L°° drift (for which (—A) 1 / 2 is critical). 



Proof of Theorem [775] The proof consists of applying Theorem 3.8 to the drift-diffusion equation 
( |7.1[ ). Since we wish to consider drift in L°°(IR 2 ) it is natural to take the same drift as in Section]!} 
but which we smoothen as done in Section [3] The main difficulty arrises in proving condition 
p.22[) in Theorem |3.8[ The issue is that by considering C, we have lost the homogeneity of the 



associated kernel. Most of the analysis is devoted to finding a constant H(a) such that p.22[ ) holds 
for all a > 0. 

Uo ;P (x), which is defined by letting a = 
norm independent of p (which 



Since we work with bounded drift, we set u p (x) 



(and r a = 1/4) in ( |3.4| ). In particular, u is divergence-free, has L c 
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we normalize to be less than 1), is globally smooth and vanishes in a ball B Ep m (0), where e Ptm will 
be chosen later on depending on p and m. In addition, for x £ C p = (B 100 (0) \ B £p m (0) ) D (C U C ) 
we have the explicit formula 



u 



p.m. 



[x 1 ,x 2 ) = sgn(a;2)(0, -1] 



(7.5) 



where C is the cone centered at the origin which is tangent to Bi/ 2 (0, 1), and C' is the reflexion of 
Co about the origin. In particular 



u{0,x 2 ) = (0,-1) 



(7.6) 



for all x 2 £ (£p,m, 100). Thus w p m satisfies the conditions of Theorem 3.8 



As in Section [3j since the above drift obeys the symmetries of the problem and is smooth, we 
have that ( |7.1[ ) has a comparison principle on the upper half-plane, and solutions are odd in x 2 . 

Let r] E C°°(M) be a radially non- increasing, cutoff function supported on \x\ < 1, with r](x) > 
if \x\ < 1, and rj(x) = 1 if \x\ < 1/2. We set 



(f){x 1 ,x 2 ) = r](A\(x 1 ,x 2 - 1)|) - ■q(A\(x 1 ,x 2 + 1)|). 



(7.7) 



Without loss of generality, we may assume that 9 > 4>{x) in the upper half-plane {x 2 > 0}. As 
outlined in Theorem 3.8 we define z(t) as the solution of 



i(*)=U2(0, *(*)) = OW), z(0) = l, 



which is given explicitly as 



for all i: G [0, t p>m \. Here t 



z(t) = l-t 



(7.8) 



p,m 



1 — £ P ,m is the time it takes z(t) to reach the value e P:Tn (we 
need to work in this time interval in order to use the explicit formula ( |7.6[ )). Since the vector 
field considered here is the same as for a = in Section [3} condition ( 3.23[ ) of Theorem |3.8 
automatically holds by p.!3[ ) in the proof of Lemma |3lT 
For any a > 0, we define 



[x 



(x/a), where is defined in (|7.7[ ). We now need to quantify 
the effect of C on a . Verifying that condition ( |3.22| ) of Theorem 3.8 holds for some constant H(a) 
is more delicate than in Section |3j so we state this in a lemma, which we shall prove at the end of 
this section. 



Lemma 7.1 (Effect of dissipation). For any a e [0, 1] we have that 

£0aO) < H(a)<f) a (x) 
holds in the upper half plane {x 2 > 0} where 



H(a) 



Co 



rm(r) dr + 



1 m(r) 



dr 



(7.9) 



(7.10) 



for some positive constant c which may depend on r\ and m, but not on a. 

We conclude the proof of Theorem |1.5[ and then return to prove Lemma 7.1 It is left to verify 
condition (|3.24|) in Theorem |3.8[ For this purpose we first prove that 



H(z(t))dt < C < oo, 



(7.11) 
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for some constant C that is independent of e p>m (and hence independent of p). Note that once 
( |7.11[ ) is proven, we have that at time T = t Pjm — 1 — c 



- p,m 



p{2z{T)) = p(2e p>m ) < exp(-C7) < exp - / H(z(t)) dt 



by choosing e„ m sufficiently small, thereby proving condition (3.24) in Theorem 3.8 



We now prove ( |7.11| ). By ( |7.8| ) and ( |7.10[ ), since c does not depend on £ Ptm , and since t p>m < 1, 
it is sufficient to estimate 



Z (ty 
i 



z(t) 



rm(r) dr + 



mir 



dr dt 



*(*) 



< 



o 



(I-*) 1 



l-t 



rm(r) dr + 



1 m{r) 



l-t 



dr dt. 



First, using Fubini we have 

'"' rL m{r) 



dr dt 



o J l-t 



1 m(r) 



ldt dr 



l-r 



m(r) dr < oo. 



Second, cf. 7.3 we use that sm(s) is non-decreasing on (0, 1), and obtain 



(l-t)* 



l-t 



rm(r) dr dt < 



[ - — ] —r [ (l-t)m(l-t)drdt= [ m(l - 1) dt < oo 
Jo I 1 - t) Jo Jo 



which concludes the proof of ( |V. 1 1 [ ). 



This concludes the proof of Theorem 1 .5 modulo Lemma 7.1 which we prove next. 



□ 



Proof of Lemma 7.1 The proof is similar to that of Lemma 3.4 but several difficulties arise be- 
cause we lost the homogeneity of the kernel of C. 

Recall that supp<p a H {%2 > 0} = B a u(0,a) =: Vt a . We prove prove that on dCl a , we have 
that £(f) a < —5 a < for some 5 a > 0. The reason this holds is that point on dVl a are local 
minima of O , and the positive contribution from the lower half plane is dominated by the negative 
contribution from the upper half plane. Let x E dQ a . By the definition of <\> in ( |7.7| ) we have that 
4> a (x) =0 and therefore 



C(f) a (x) 



m(\x-y\) 

My)-, 

\x — y\ z 
M(yi,V2 - a)\\ m(\x-y\) 



\x - y\ 



dy + 



M(yi,V2 + a)\\ m(\x-y\) 

^dy 



\x - y\ 



n 



M(yi,V2 -a)\ \ ( m(\x-y\) _ m(\x - Ry\) \ 



\x - y\ 



\x — Ry\ 2 J 



dy 



(7.12) 



where for y = y 2 ) we have denoted Ry = (y x , —y 2 ). For x E dVt a and y E f2 a , we have that 
\x — y\ < \x — Ry\, and hence due to the monotonicity of m, the integrand in ( |7.12[ ) is positive. 
Coupled with the fact that 77 > 0, this already shows Ccj) a (x) < 0. Note that for x E <9-B a / 4 (0, a) 
and y E -B a /4.(0, a), by the triangle inequality we have \x — Ry\ > 2\x — y\. This can be seen by 
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drawing a picture. Therefore, since m is decreasing we obtain from ( |7.12| ) that 

'4\(yi,y 2 - a)\\ f m(x - Ry) m(x - Ry) 



C(f) a (x) < 



V 



B o/8 (o,a) v a / V \x-y 
3 f m(3a) , m(3a) 



< 



4|x — y\ 



dy 



4 



S o/8 (0,a) \ x V\ 



~5a 



(743) 



and (17. 91) holds for x G <9f2 Q 



Similarly, since |x — y\ < \x — Ry\, whenever both x and y are in the upper half-space, on the 
set {x 2 > 0} n fi„ we have £0 a < 0, so that (L9\ trivially holds. 

By smoothness of 77, there exists p a G (0, a/4) such that on the annulus B a u(0, a)\B a / i _ Pa (Q, a) 
we have C<p a < 0. We have to estimate p a from below, as this will be needed later on. Let 
x £ -B a / 4 (0, a) \ -B a /4_ Pa (0, a) for some p a > 0. From the mean value theorem, setting 



p a = mm 



S„ 



\\VC(/) a \\ L <x> ' 

ensures that C<p a (x) < —S a /2. In addition, we have 

|V£0 o (a;)| < / |V0 fl (x) - V0 a (x + y)| ^ 



(7.14) 



dy 



< c\\V<j) a \ 



C 1 



\v\< 



™>(y) 

1 \y\ 



dy + c||V0a|U- 



\y\>l \y 



m(y) 

— 2~ d y 



< 



'll'/llr-- A 1 , s , , ^Nb rmyr) 

m(r)dr-\ / dr 



a 







a 



r 



(7.15) 



in view of (7.2). It follows from (7.14) and (7.15) that 



since m is monotone decreasing and a < 1. Here C3 < 1/2 is independent of a. 



(7.16) 



Lastly, if x G _B a /4_ Pa (0, a), then since 77 is radially non-increasing, by ( |7.16| ) we have 

/ 4|(xi,x 2 -a) \\ ^ 4p a ' 



a (x) = 77 



>77 1- 



> > 



(7.17) 



where c v is independent of a. Therefore, to ensure that C(j) a (x) < H(a)4> a (x), we just need to 
verify 



H(a) > 



L°° 



(7.18) 



Similarly to ( |7.15[ ) we may bound (but this time we split the integral domains at \y\ = a not at 
\y\ = 1, and we exploit the P.V. in the definition of C to write the nonlocal operator in terms of 
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double-differences) 



\C(p a {x)\<- / \2<p a {x)-<p a {x + y)-<p a {x-y)\-^dy 
2 J®? \y\ 2 

^ c\\r]\\ C 2 f a f 1 m(r) , f°° m(r) . 
< — / rmir) dr + c / dr + c / dr 



° 2 -h Ja r 



< '-a I -4 f rm(r)dr+ [ dr + 1 ) (7.19) 



for some c > independent of a. Combining (7.18H7.19) completes the proof of the lemma. □ 
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